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We find that if two superconducting isiands of different number parity are finked by a tunnei junction 
the unpaired eiectron in the odd isiand has a tendency to tunnef into the even isiand. This process 
ieads to ffuctuations in time of the number parity of each isiand, giving rise to a random teiegraph 
noise spectrum with a characteristic frequency that has an unusuai temperature dependence. This 
new phenomenon shouid be observabie in a Cooper-pair pump and simiiar singie-eiectron tunneiing 
devices. 



PACS numbers: 74.20.Fg, 74.40.+k, 74.50+r 

In this paper we study the problem of tunneling be- 
tween number restricted superconducting islands. Small 
capacitatively isolated superconducting islands have re- 
cently attracted much experimental [1] and theoretical 
[2,3] attention. The basic question addressed in vari- 
ous ways by these studies is: how are the properties of 
a superconducting sample affected by the evenness or 
oddness of the number of electrons it contains? Since 
the Bardeen-Cooper-Schrieffer (BCS) [4] ground state is 
made from a coherent superposition of Cooper pairs there 
will always be an unpaired particle present if the total 
number is odd. In the experiments cited the effects of 
number parity are revealed by using the Coulomb energy 
to control the addition or removal of a single electron. 

When an unpaired particle is present, the energy dif- 
ference between an odd and an even sample at absolute 
zero will always be larger by an amount A than it is in the 
normal state: E(2N +1) — E(2N) — /j, = A no matter how 
big N is (here A and /j, are the superconducting energy 
gap and chemical potential, respectively). At low, but fi- 
nite temperatures, the free energy difference between the 
two samples is equal to S^odd-even = A — fc^TTn N e ff. 
The entropy term —ksTlnNeff arises because the un- 
paired particle can sample the available states (approx- 
imately N e ff in number). It is only a logarithmic func- 
tion of N since N eff ~ 2A(0)VV2ttTA oc N [1-3]. 
[Here N(0) is the single particle density of states at 
the Fermi level, and V is the volume of the island.] At 
T* = A/ In N e ff ~ 200 - 300mA' the free energy differ- 
ence practically vanishes (for typical island parameters 
see Ref. [1]). As a result parity effects are robust against 
size; they have been observed [1] in the equilibrium and 
transport properties of islands with particle numbers as 
large as roughly N ~ 10 9 . The constraint on N comes 
from the Coulomb blockade condition of thermal fluctu- 
ations k B T < e 2 /2C. 

Now consider two islands between which electrons can 
tunnel. The particle number and parity of each individ- 
ual island is no longer fixed, although the total number 
and parity of the system is conserved. The situation 
greatly simplifies at low temperatures, when very few 
quasiparticles are present. Cooper pair tunneling, which 



does not affect the parity of the island, corresponds to a 
twice as large charge fluctuation and can be suppressed 
by the Coulomb blockade even when the charging energy 
of the system is tuned by gate voltages to permit single 
particle tunneling. The most interesting case is when ini- 
tially one of the islands is in the even, and the other one in 
the odd parity state. When T < T* , all the particles are 
paired up on the even island, and on the average there 
is only one quasiparticle present on the odd island. If 
the tunnel resistance R T > R Q = h/e 2 [5], the unpaired 
electron is well localized on the odd island. Nevertheless, 
there is a nonzero transition rate for tunneling of a sin- 
gle quasiparticle, always from the odd island to the even 
island. This leads to a random telegraph noise type fluc- 
tuation in time of the number parity, with a Lorentzian 
noise spectrum [6]. 

An important remark should be made at this stage. 
This parity oscillation in this regime is not a coherent 
quantum oscillation: the state of the unpaired particle is 
not given by a phase coherent superposition of its local- 
ized wave functions on the left and right islands. Because 
there are many available states for the particle on each 
island the transition rate must be evaluated by Fermi's 
Golden Rule and no phase coherence can be preserved 
under these circumstances. This incoherent tunneling 
process, however, as we will show later on, has some very 
peculiar properties: it sets in below T* and the tunneling 
rate either has a maximum or saturates as the temper- 
ature is lowered, depending on how strong the coupling 
is between the system and the environment. Such a be- 
haviour is quite unusual and cannot be easily sorted into 
the any of the usual classes, like that of thermally acti- 
vated processes [7]. 

We will suggest experimental setups in which this spec- 
trum could be measured. Our analysis gives new stability 
regions for the Cooper pair pump [8] and other devices af- 
fected by parity effects and implies modified parameters 
for optimal operation. The coupling between the tunnel- 
ing electron and the collective modes of the electromag- 
netic environment turns out to be very important. Our 
calculations show that when this coupling is too strong 
the unpaired electron becomes trapped on one island. 
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Our model consists of two superconductors of different 
parity, even and odd, linked by a tunnel junction. For 
the sake of simplicity, we take these superconductors to 
be otherwise identical: the same volume V, energy gap 
A etc. We want to describe the tunneling events tak- 
ing place in this junction. Since electrons have charge, 
a tunneling event perturbs the equilibrium charge distri- 
bution, and excites environmental modes: oscillations of 
the electromagnetic field in the circuit (or electromag- 
netic environment) connected to the junction. Tunneling 
in the presence of environment is assumed to be ade- 
quately described by the Hamiltonian [9,10]: 
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The first term describes the tunneling with amplitude 
Tkq of an electron in state ka from the left into qa on 
the right. The second terms corresponds to the the time- 
reversed process. We assume zero magnetic field and 
make the arbitrary choice that the superconductor on 
the left is odd and the one on the right is even in parti- 
cle number. A e is a charge displacement operator which 
satisfies A e QA\ = Q — e. The quantity of interest is 
the transition rate of the system from this initial state 
into one in which the parity of the superconductors is 
switched. The total rate is the difference of quasiparti- 
cle tunneling rate from left to right, r _> (V), and from 
right to left, r* _ (V). At low temperatures and zero bias 
voltage this rate, in contrast with the usual situation, 
is nonzero and is directed from the odd superconductor 
towards the even one (in this case from left to right): 
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where R^ 1 = Ait(e 2 / h)\T\ 2 , and Rt is the tunneling re- 
sistance. It is assumed that |X)% g | 2 is a weak function 
of momentum and can be replaced by its average \T\ 2 
around the Fermi energy. The quasiparticle density of 
states is given by the BCS result [4] 
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and zero otherwise. The function P(E) is the Fourier 
transform of (At(t)A(O)} [10]. In physical language, 
P(E)dE gives the probability that the electron exchanges 
an energy in the range E, E + dE with the environment 
during the tunneling process. 

A similar process occurs in the ordinary quasiparticle 
tunneling [11] when thermally excited quasiparticles tun- 
nel between macroscopic superconductors giving rise to 



a structure in the I-V curve at V = ±|Ai — A2I. In this 
case there are no thermally excited quasiparticles but 
there is one available purely due to the parity constraint. 
The difference in the distribution functions f — f e which 
enters the above expression, is [3] 
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At low temperatures the difference in quasiparticle oc- 
cupation numbers becomes f — f e = exp[— f3(Ek — 
A)]/2N e ff whereas close to T* the denominator in (4) 
becomes very large and the parity difference becomes 
negligible. 

While it is not possible to compute P(E) in a general 
case, there are two important and physically intuitive 
limits to consider. For low environmental impedance 
|Z(w)| <C h/e 2 the dominant process is elastic tunnel- 
ing: P(E — E') = 8(E — E'). The environmental charge 
configuration remains the same, similar to the recoilless 
Mossbauer transition. If, however, |Z(w)| h/e 2 , at low 
temperatures T <C E c , P(E) becomes 



P(E) = 8(E-E C ), 



(5) 



where E c is a charging energy scale (note that an inter- 
polation between the two limits is possible formally turn- 
ing E c to zero). In the latter limit for low temperatures 
T < T* < E c we obtain 



r= Aexp(-^ e)(l + e . 



f)Aa 



2e 2 R T N eJJ 



+00 



(x + a)xe 



-f)A(x-l) 



V(* 2 -l)[(z + a)2-l] ! 



(6) 



where a = E c /A. When a = the integral diverges 
logarithmically, as does the ordinary quasiparticle cur- 
rent between two macroscopic superconductors [11] due 
to the overlapping BCS density of states. It is quite well 
understood for a long time [12] that in a real metal the 
divergence is naturally removed by gap anisotropy and 
nonzero quasiparticle lifetime. For nonzero a the rate is 
further suppressed by the overal factor exp(— /3E C ). The 
odd electron is "locked" by the electromagnetic environ- 
ment into one island: a tunneling event in the presence 
of a high impedance environment would result in the 
creation of a large number of low energy environmental 
modes so that the initial environmental state is practi- 
cally orthogonal to the required final state. As a con- 
sequence, the transition amplitude is decreased by this 
"orthogonality catastrophe" [13]. 

The dependence of the fluctuation rate on other pa- 
rameters is clear from the above result and easy to inter- 
pret. The process is slowed down by increasing tunnel 
resistance, which corresponds to decreasing the trans- 
parency of the junction. The time spent by the odd 
particle in one superconductor, T _1 is proportional to 
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the volume of the superconductor via N e j j , as it should 
be. In order for the process to be observable, the tem- 
perature has to be low enough, T < T* , so that parity 
effects are important. Below this temperature the rate is 
proportional to the energy gap. 

The complete temperature dependence of T(T) is 
shown in figure 1, calculated from eqs. (2), (4) and (5). 
First note that the rate becomes nonzero only below T* 
and has a maximum before it becomes exponentially sup- 
pressed according to eq.(6) due to the coupling to the en- 
vironment. The parameters were chosen from reference 
[8] for later convenience, when we discuss these ideas in 
the context of the Cooper pair pump. If the coupling is 
weak (see inset), the rate levels off to a very high value 
limited by gap anisotropy or impurity scattering. These 
effects smear out the square-root singularity in the BCS 
density of states. We treat this limit by formally taking 
the parameter a to a very small but nonzero value. Such 
temperature dependence is quite unusual for incoherent 
quantum tunneling process [14] where the exponential 
dependence (corresponding to thermally activated hop- 
ping) at higher temperatures is taken over by power law 
dependence T oc J 12 " -1 at low temperature. 

Let us now turn to the discussion of the single-Cooper- 
pair pump [8]. In our opinion, it has all the necessary 
ingredients to exhibit parity fluctuations. It is the sim- 
plest device that has two distinct islands (see Fig. 2), 
and can be charged in a locally stable manner - by ap- 
propriately tuning the gate voltages Ui and U2 - to a 
configuration characterized by (rii,?^)- (Here n\ and 
«2 are the number of excess electrons on each island). 
We have calculated the regions of stable charging in the 
(Ui, U'z) plane. Our result is different from the previous 
ones [8] in that we now take into account the increase 
in the free energy of those configurations which unpaired 
electrons are present. In each such region the system has 
to be stable against tunneling one each junction i to right 
(left) i+, (i-) (here i = {1, 2, 3}) 

AF {± = AE i£ + A[p(ni) + p(n 2 )].T e/o > 0, (7) 

That is, the change in the total free energy of the system 
AF l:iz equals to the sum of changes in the charging en- 
ergy AE l:iz and parity dependent energy, 6[p(rii)+p(n2)], 
where p(n) = 1 if n is odd and p(n) = if n is even. As 
the temperature is lowered below T* the regions with odd 
island charge configuration will shrink, especially those 
with both n\ and odd. The result of the calculation 
is shown in Fig. 3 for a particular value of the parameter 
d = T e j j{e 2 /3Co) = 1/2 when the even-odd boundaries 
(dotted lines) have a maximal length for nonzero T e / . 
(Co is the capacitance of the tunnel junctions, assumed 
to be equal). At this stage the all-odd configurations, 
like (1, 1), disappear from the ground state configura- 
tion. We obtain a very simple condition necessary for 2e 
charge quantization: A > e 2 /3Co- This is desirable for 



the successful operation of the Cooper-pair pump since 
then only the all-even regions survive the low temper- 
ature limit. The condition seems not to be fulfilled in 
ref. [8] where even at the lowest operating temperatures 
(50mA') d ~ 0.9. 

Assuming that both islands have an even number of 
electrons in the configuration (0,0), we are interested in 
the case of n\ odd and even (or vice versa). The sys- 
tem must to be tuned via JJ\ and U2 at the lines separat- 
ing the stability regions (2n+l, 2n) from (2n, 2n+l) (dot- 
ted lines in Fig 3 ). The circle shows a typical load curve 
around which the Cooper pair pump is driven. Clearly, 
the system is unstable against parity fluctuations. An 
error in the current occurs if such an event forces the 
system to leave the circle. Previous studies [8] of the 
Cooper-pair pump reported difficulties connected to the 
presence of quasiparticles at very low temperatures. Par- 
ity effects might be at the origin of these difficulties, since 
an island with odd number of particles contains at least 
one quasiparticle, no matter how low the temperature is. 

By coupling the islands to single-electron electrome- 
ters, the process should be detectable by measuring the 
island charges as a function of time and extracting the 
noise spectrum. Above T* there should be no effect, and 
the values of n\ and should evolve relatively smoothly 
from one value to the other one, as the system is driven 
accross such boundary. Assuming, as is reasonable, that 
there is no correlation between the parity flip events and 
that the probability for such a transition in the time 
interval t,t + dt is T(T)dt, the parity fluctuation pro- 
cess will generate telegraph noise [6] with its character- 
istic Lorentzian noise spectrum S(ui) = \Y{lu 2 + AT 2 )' 1 . 
Focusing on the middle junction which controls parity 
fluctuation between the islands, we approximate the cir- 
cuit by replacing the outer junctions with capacitors [9]. 
Then the charging energy scale used in eq. (5) becomes 
E c = l/3(e 2 /2C ) and our result (6) for T(T) can be 
directly used. Furthermore, we assume that the time be- 
tween each tunneling event is larger than the characteris- 
tic relaxation time, so that it is justifiable to use the equi- 
librium environmental states in calculating the tunneling 
rates. We have not considered possible quasiparticle tun- 
neling events on one of the outer junctions since they can 
be neglected in a system with significant Coulomb block- 
ade. Higher order processes, like co-tunneling, as well as 
quantum fluctuations of the island charge are only im- 
portant when the tunnel resistance is comparable to the 
resistance quantum. 
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FIG. 1. The temperature dependence of the parity 
fluctuation rate. For high environmental impedance 
|^(a;| ^ h/e 2 . The inset shows the case when 
|^(a;| <C h/e 2 . (parameters from Geerlings et al.: 
T* ~ 220mA", N eff ~ 10 4 , R T = 85A;Q, A = 2. 31 A", a = 
E c /A = 0.42, see eq(7) and text). 




Figure 2. 

FIG. 2. Schematic circuit diagram of the single Cooper pair 
pump, with excess electrons n\ and «2 on islands and delim- 
ited by dashed lines. Gate voltage sources U\ and U2 are 
coupled to the islands via two capacitors. The parity fluctua- 
tion is mediated by the middle junction, and is monitored by 
the electrometers E\ and E2. In this paper we take V = 0. 




Figure 3. 

FIG. 3. Stability domains for the ground state charge con- 
figuration of the single Cooper-pair pump, in the plane of gate 
voltages {1/1,1/2). Numbers (rai,W2) correspond to the num- 
ber of excess electrons on each island. Thin dashed lines give 
boundaries for stable single electron charging, whereas thick 
solid lines give the Cooper-pair stability regions. The dot- 
ted lines mark the boundary where parity fluctuation should 
be observed. The circle is a typical cycle load curve used to 
operate the Cooper-pair pump 
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